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Abstract — On a two-level quantum system driven by an 
external field, we consider the population transfer problem from 
the first to the second level, minimizing the time of transfer, 
with bounded field amplitude. On the Bloch sphere (i.e. after 
a suitable Hopf projection), this problem can be attacked with 
techniques of optimal syntheses on 2-D manifolds. 

Let {—E,E) be the two energy levels, and \fl{t)\ < M the 
bound on the field amplitude. For each values of E and M, we 
provide the explicit expression of the time optimal trajectory 
steering the state one to the state two in terms of a parameter 
that should be computed numerically. 

For M « E, every time optimal trajectory is periodic (and 
in particular bang-bang) with frequency of the order of the 
resonance frequency lor — 2E. 

On the other side, for M > E the time optimal trajectory 
steering the state one to the state two is bang-bang with exactly 
one switching. Fixed E we also prove that for M ^ oo the time 
needed to reach the state two tends to zero. 

Finally we compare these results with some known results of 
Khaneja, Brockett and Glaser and with those obtained in the 
Rotating Wave Approximation. 

Keywords: Control of Quantum Systems, Optimal Syn- 
thesis on the Bloch Sphere, Minimum Time 

I. Introduction 

In this paper we apply techniques of optimal synthesis on 
2-D manifolds to the population transfer problem in a two- 
level quantum system (e.g. a spin 1/2 particle) driven by an 
external field (e.g. a magnetic field along a fixed axis). Two- 
level systems are the simplest quantum mechanical models 
interesting for applications (see for instance [3], [10]). The 
dynamics is governed by the time dependent Schrodinger 
equation (in a system of units such that fi = 1): 
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where E, is a real number {±E represent the energy levels 
of the system). The control ri(.), that we assume to be a 
real function, different from zero only in a fixed interval, 
represents the external pulsed field. In the following we call 
drift term , the Hamiltonian with no external fields (i.e., the 
term diag{~E, E)). 

The aim is to induce a transition from the first level (i.e., 
iV'iP = 1) to the second level (i.e., |V'2p — 1), minimizing 
the time of transfer, with bounded field amplitude: 



where T is the time of the transition and M is a positive real 
constant representing the maximum amplitude available. 

Remark 1: This problem was studied also in [12], but with 
quadratic cost fl{t)^ dt and with no bound on the control. 
In this case, optimal solutions can be expressed in terms of 
Elliptic functions. 

It is a standard fact to eliminate an irrelevant global factor 
of phase by projecting the system on a two dimensional real 
sphere S*^ (called the Bloch Sphere) by means of an Hopf 
map. In this way the Schrodinger equation (1), (2) becomes 
the single input affine system (after setting u{t) ~ il{t)/AI): 
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Fs{y) +uGsiy), where: 
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2/= (yi,y2,2/3) GK^ = 
\u\ < 1, 

-2/2 

Fs{y) ■■= fccos(a) I 2/1 
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a = arctan ( — j e ]0, 7r/2[, 
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and the constant k is given by k 



2E/cos{a) 



\n{t)\ < M, for every / e [0, T], 



(3) 



Normalizations In the following, to simplify the notations, 
we normalize fc = 1. This normalization corresponds to a 
reparametrization of the time. More precisely, if T is the 
minimum time to steer the state one to the state two for the 
system with fc = 1, the corresponding minimum time for 
the original system is simply T/V AI^ + E^. We come back 
to the original value of fc only in Section III-C. 

The vector fields Fs{y) and Gs{y) describe rotations 
respectively around the axes ya and j/i. Now the state one is 
represented by the point 2/3 = 1 (called the north pole) and 
the state two by the point 2/3 = —1 (called the south pole). 
The optimal control problem is then to connect the north 
pole to the south pole in minimum time. As usual we assume 
the control u{.) to be a measurable function satisfying (6) 
almost everywhere. The corresponding trajectory is a 
Lipshitz continuous function y{.) satisfying (4) almost 
everywhere. 

The most important and powerful tool to study optimal 
trajectories is the well known Pontryagin Maximum Principle 
(in the following PMP, see for instance [2], [13]). It is a first 



order necessary condition for optimality and generalizes the 
Weierstrafi conditions of Calculus of Variations to problems 
with non-holonomic constraints. For each optimal trajectory, 
the PMP provides a lift to the cotangent bundle that is a 
solution to a suitable pseudo-Hamiltonian system. 
Even if the PMP is powerful, giving a complete solution to 
an optimization problem remains extremely difficult. First, 
one is faced with the problem of integrating a Hamiltonian 
system. Second, one should manage with "non Hamiltonian 
solutions" of the PMP, the so called abnormal extremals . 
Finally, even if one is able to find all the solutions of the 
PMP it remains the problem of selecting, among them, the 
optimal trajectories. For these reasons, usually, one can hope 
to find a complete solution of an optimal control problem for 
very special costs, dynamics and in low dimension only. 

Two dimensional minimum time problems in control affine 
form, (like the problem (4)-(9)) are nice cases for which the 
analysis can be pushed much further, thanks to the theory 
developed in [7] (see also the references therein). In this 
paper we take advantage of that theory to restrict the set of 
candidate optimal trajectories. The optimal trajectories are 
then identified, by requiring that they respect certain crucial 
symmetries of the system. More precisely, for M > E the 
time optimal trajectories steering the state one to the state 
two are bang-bang with exactly one switching, and we give 
the exact expressions of the corresponding optimal controls. 
In particular, fixed E we see that for M oo the time 
needed to reach the state two tends to zero. 
On the other side, for M << E, every time optimal 
trajectory is periodic (and in particular bang-bang) with 
frequency of the order of the resonance frequency ojji ~ 2E, 
and can be selected among a finite set of trajectories which 
corresponds to solutions of suitable equations. 

Remark 2: If we were describing either a spin 1/2 particle 
driven by two magnetic fields (one along the x axis and 
one along the y axis) or a two-level molecula driven by an 
external field in the Rotating Wave Approximation (RWA 
for short, see for instance [3]), then our Hamiltonian would 
contain complex controls: 
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where (*) indicates the complex conjugation involution. In 
this case the minimum time problem with bounded controls 
(i.e., |ri(t)| < M) is easier, since it is possible to eliminate 
the drift term by a unitary change of coordinates and a change 
of controls (interaction picture). This problem has been 
studied in [8], [12]. The simplest time optimal trajectory, 
steering the system from the state one to the state two, 
corresponds to controls in resonance with the energy gap 
2E, and with maximal amplitude i.e. 



(11) 



The quantity ojr = 2E is called the resonance frequency. 
In this case, the time Tfjw-^ of transfer is proportional to 
the inverse of the laser amplitude. More precisely T^wyi — 
7r/(2M), see for instance [8]. In Section (III-C) the minimum 



time of transfer for the Hamiltonians (2) and (10) are 
compared. 

In Section II, we recall some basic properties of optimal 
trajectories for the system (4)-(9), that were already obtained 
in [5]. In Section III we state our main results and in III-C we 
compare these results with some known results of Khaneja, 
Brockett and Glaser and with those obtained in the Rotating 
Wave Approximation. In Section IV we give an idea of the 
techniques we used. 

II. Known Results 

Recall that we have normalized k = 1. Note that, since 
the system (4)-(9) is Lie bracket generated on a compact 
manifold and the set of velocities is compact and convex 
then, for each pair of points p and q belonging to 5^, there 
exists a time optimal trajectory joining p to q. 

Although, with different purposes. The minimum time 
problem for the control system (4), (9), although with dif- 
ferent purposes, has been partially studied in [5] (see also 
[6]). 

In particular in [5] it was proved that every time optimal 
trajectory is a finite concatenation of bang arcs (i.e., cor- 
responding to control a.e. constantly equal to +1 or —1) 
and singular arcs (i.e., corresponding to singularities of the 
End point mapping, see for instance [4], that in our case 
correspond to controls a.e. vanishing) with some special 
structure. More precisely: 

Definition 1: A control u : [a, b] [—1, 1] is said to be 
bang-bang if u{t) £ {—1, 1} a.e. in [a, b]. Moreover, ifu{t) E 
{ — 1, 1} andu{t) is constant for almost every t S [a, b], thenu 
is called a bang control. If u{t) = {) for almost every t G [a, 6], 
then u is called a singular control. 

A switching time ofu is a time t G [a, b] such that, for every 
e > 0, u is not bang on {t ~ e,t + e) O [a, b]. A control with 
a finite number of switchings is called regular bang-bang. A 
trajectory of the control system (4)-(9) is a bang trajectory, 
singular trajectory, bang-bang trajectory, regular bang-bang 
trajectory respectively, if it corresponds to a bang control, 
singular control, bang-bang control, regular bang-bang control 
respectively. 

In the sequel, we use the following convention. The letter 
B refers to a bang arc and the letter S refers to a singular 
arc. A concatenation of bang and singular arcs is labeled by 
the corresponding letter sequence, written in order from left 
to right. Sometimes, we will use a subscript to indicate the 
time duration of an arc, so that we use Bt to refer to a bang 
arc defined on an interval of length t and, similarly, St for 
a singular arc defined on an interval of length t. 

Using the PMP and the Theory developed in [7], in [5] 
(see also [1]) it was proved the following: 

Proposition 1: Consider the control system (4)-(9). Then: 

A. every time optimal trajectory is a finite concatena- 
tion of bang and singular arc; 

B. if a time optimal trajectory contains a bang arc Bt, 
then t <2tt; 
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Fig. 1. Graph of v(.) when a = 7r/6 



C. if a time optimal trajectory contains a singular arc, 
then it is of the type BtSgBf, with s < ^^J^^^ , 
t, t' > 0. Moreover the support of singular arcs lies 
on the set (called equator) = 0. 

D. if a time optimal trajectory is bang-bang, then the 
time duration T along an interior bang arc is the 
same for all interior bang arcs and verifies n <T < 
27r. 

From C. it follows that the first and the last arc on optimal 
trajectory connecting the north with the south pole are not 
singular. The following proposition (see [5] for the proof) 
gives more details on the optimal trajectories starting at the 
north pole in the case a < tt/A. 

Proposition 2: Consider the control system (4)-(9), and 
assume a < tt/4. Then the optimal trajectories starting from 
the north pole are of the form Bg^ By(^Si) • ■ • By(^Si)Bsf , where 
Si e [0,7r], sf G [0,w(sO] and 



v{si) = TT + 2 arctan 



sin(Si) 



(12) 



^cos(si) + cot^(a) ^ 
Note that the function v{.) is such that v{0) = v{tt) = it 
and moreover it is increasing on the interval [0,t\ and 
decreasing on [t, tt], where i ~ arccos(— tan^(a)), 
moreover if a is small the maximum of v{.) is 
v{t) = 2 arccos(— tan^(Q;)) ~ tt + 2a^ (see Figure 1). 



III. Main Results 

A. The a > tt/4: Case 

In the case a > tt/A, there are exactly four optimal 
trajectories steering the state one to the state two. They are 
easily described by the following: 

Proposition 3: Consider the control system (4)-(9), and 
assume a > tt/A. Then the optimal trajectories steering 
the north pole to the south pole are bang-bang with only 
one switching. More precisely they are the four trajectories 
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Fig. 2. Case a > 7r/4: optimal trajectories 

corresponding to the four controls 

1, te [0,sa] 
-1, te\sA,T] 

,(3) ^ J -1, t e [0, sa] 

1, t^]sA.T] 



,(2) 



,(4) 



1, t e [o,sb] 

-1, t£]sB,T] 

-1, i e [0,ss] 
1, t <e]sb,T] 



where: 



sa^ — arccos(cot^(Q!)), 
SB = TT + arccos(cot^(Q;)), 
T = 2tt. 



(13) 



One can easily check that the switchings described in Propo- 
sition 3 occours on the equator (ys — 0). Figure 2 represents 
the optimal controls u^^^ and u^^). 



B. The a < tt/A Case 

If a < tt/A, the situation is more complicated. From 
Proposition 2, we know that every optimal trajectory starting 
at the north pole has the form: 



Bs. B. 



(14) 



n— 1 times 



where the function v{si) is given by formula (12). (In the 
following we do not specify if the first bang corresponds to 
control +1 or —1, since, as a consequence of the symmetries 
of the problem, if u{t) is an optimal control steering the 
north pole to the south pole, —u{t) steers the north pole to 
the south pole as well.) 



Fig. 3. Graph of the function T with a = 0.13 



Fig. 4. Graph of the function Q with « = 0.13 



It remains to identify one or more values of Si , s/ and the 
corresponding number of switching n for this trajectory to 
reach the south pole. 

Next, given s e [0, tt] such that s ^ t = 
arccos(— tan^(a)) we call s'(s) the unique solution to the 
equation v(s) = u(s'(s)) with s'(s) 7^ s and we define 
s'(t) = t (see also Figure 1). Considering the symmetries 
of the problem, one can prove that if a < 7r/4, s/ is equal 
either to Si or to s'(si). This fact is described by Lemma 1 
below. 

In the following we describe how to identify candidate 
triples {si,Sf,n) for which the corresponding trajectory 
steers the north pole to the south pole in minimum time. 
There are two kind of candidate optimal trajectories. 

• Sf = s'{si), called TYPE- 1 -candidate optimal trajecto- 
ries 

• Sf ~ Si called TYPE-2-candidate optimal trajectories 



Define the following functions: 

v{s) 



9{s) ~ 2arccos|^sin 



cos(2q;) — cos — — 



v{s) 



(15) 



(3{s) = 2 arccos(sin(a) cos(q!)(1 — cos(s))) (16) 

Proposition 4: (TYPE-l-trajectories) Let a < n/i and 

s £ [0,7r]. Fixed a, the following equation for the couple 

(s,n); 



27r 



(17) 



has either two or zero solutions. More precisely if (s,n) 
is a solution to equation (17), then (s'(s),n) is the sec- 
ond one, and the trajectories Bg B^i^gy ■ ■ B^^^g^ B^'i^s) ^nd 

" V ' 

n-l 

Bs'{s) Bv(^s) • • • Su(s) Bg are TYPE- 1 -candidate optimal tra- 
^ ^ ■' 

jectories. 

In figure 3 the graph of the function (17) is drawn for a 
particular value of a, namely a ~ 0.13. 

Proposition 5: (TYPE-2-trajectories) Let a < n/i and 

s e [0,7r]. Fixed a, the following equation for the couple 
(s,n): 



y(s) „, . +l = n, 



has exactly two solutions. More precisely these solutions 
have the form {si,n), {32,1^ + 1) and the trajectories 

Bsi By(si) ■ ■ ■ By(si) Bsi and Bg^ -Su(s2) ' ' ' ^"(52) Bs2 ^re 

" V ' " V ' 

n— 1 n 

TYPE-2-candidate optimal trajectories. 

In figure 4 the graph of the function (18) is drawn for 
a = 0.13. 

Propositions 4 and 5 select a set of (possibly coinciding) 4 
or 8 candidates optimal trajectories (half of them starting with 
control +1 and the other half with control —1) corresponding 
to triples {si,Sf,n) that can be easily computed numerically. 

Then the optimal trajectories can be easily selected. Notice 
that there are at least two optimal trajectories steering the 
north to the south pole (one starting with control +1 and the 
other with control —1). 

In the particular case in which TT/{2a) is an integer number 
n one can see that TYPE-1 candidate optimal trajectories 
coincide with some of TYPE-2 candidate optimal trajecto- 
ries. They are of the type Bt^ Bt^-.-Bt^ Bt^ or of the type 

n-2 

Bs S„(s)...B^,(s) Bs for some s g]0, tt[. 



— -me [0,1[ 



denotes the integer part. One can prove the follow- 



(18) 



Otherwise if 7r/(2a) is not an integer number, define: 

r ^ 1 
Za 

where 
ing: 

Proposition 6: There exists f(m) G]0, 1[ such that: 

• if r G [0, f{m)] then equation (17) admits exactly two 
solutions that are both optimal, while TYPE-2 candidate 
optimal trajectories are not. 

• if r E]f{m), 1[ then equation (17) does not admit any 
solution. 

The claims on existence of solutions of the previous 
propositions come from the fact that J-{Q) = ^{tt) = ^ 
and the only minimum point of T occurs at 

s = TT — arccos(tan^(a)). 

It turns out that the image of is a small interval whose 
length is of order and therefore equation (17) has a 
solution only if a is close enough to for some integer 
number n. 

On the other hand it is possible to estimate the derivative 
of Q with respect to s showing that it is negative in the 



open interval ]0,7r[. Therefore, since G{0) = ^ + 1 ™d 
C/(7r) = ^ 1. equation (18) has always two solutions (if 
^ is an integer number then the trajectories corresponding 
to the solutions Si — and = tt coincide). 

Using the previous analysis one can easily prove the 
following: 

Proposition 7: If N is the number of switchings of an 
optimal trajectory joining the north to the south pole, then 



1. 



TT TT 

l< N < — 

2a ~ 2a 
Using these inequalities and the fact that the function 

2s + 1^ 11(5) is increasing on [0, tt], we can give 

a rough estimate of the time needed to reach the south pole: 

Proposition 8: The total time T of an optimal trajectory 

joining the north to the south pole satisfies the inequalities: 

2 2 

^-27T<T<^+TT. (19) 

2a 2a 

C. Comparison with results in the RWA and with [11] 

In this section we come back to the original value of k 
i.e. k = 2E/ cos(a) = 2\J Al'^ + i?^, and we compare the 
time necessary to steer the state one to the state two for our 
model and the model (in the RWA) described in Remark 2. 

For our model we have the following: 

. for a > 7r/4 then T = 27r/fc = tt/VM^ + E^; 

> for a < 7r/4 then T is estimated by 



1 . 

27r 

k\2a 



1 

k \ 2a 



On the other hand, for the model in the RWA, we have 
Trwa = 7r/(2M) (cfr Remark 2). 

Fixed i5 = 1, in Figure 5 the times T and Trwa as 
function of M are compared. Notice that although TjfwA 
is bigger than the lower estimate of T in some interval, 
we always have T^wA < T. This is due to the fact that 
the admissible velocities of our model are a subset of the 
admissible velocities of the model in the RWA. 

Notice that, fixed E = I, for M ^ we have T ^ 
TT^/ (4A/), while for M 00, we have T ^ tt/M. In other 
words: 

. for M ^ we have T - {TT/2)TjiwA, 

> for 7\jf ^ 00 we have T ^ 2TfiwA- 

Remark 3: For M << E (i.e. for a small) v{s) ^ 
tt/{2E). It follows that a time optimal trajectory connecting 
the north to the south pole (in the interval between the first 
and the last bang) is periodic with period P ^ tt/E i.e. with a 
frequency of the order of the resonance frequency ujji = 2E 
(see Figure 6). On the other side if AI > E then the time 
optimal trajectory connecting the north with the south pole 
is the concatenation of two pulses. Notice that if 7\/ >> E, 
the time of transfer is of the order of tt/M and therefore 
tends to zero as A/ ^ 00. It is interesting to compare this 
result with a result of Khaneja, Brockett and Glaser, for a 
two level system, but with no bound on controls (see [11]). 
They estimate the infimum time to reach every point of whole 
group SU{2) in n/E. 
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Fig. 5. Estimate on tlie minimum time to reacli the state two and comparison 
with the time needed in the RWA 



Indeed for our model it is possible to prove that, for AI 
00, not every point of the Bloch sphere can be reached from 
the state one in an arbitrarily small time, but this is the case 
for the state two, as we discussed above. 

IV. Sketch of the Proofs 

In this section we give an idea of the techniques we used to 
prove our results. The key point is described by the following 
Lemma which states a property of optimal trajectories as a 
consequence of the symmetries of the problem. Recall that 
we have normalized k = 1. Before stating the Lemma we 
note that it is possible to extend Proposition 2 to the case 
a > 7r/4, assuming that Si E [0, arccos(— cot^(a))[ , and 
in this case v{.) is an increasing function on its interval of 
definition. 

Lemma 1: Every optimal bang-bang trajectory, connect- 
ing the north to the south pole, with more than one switching 
is such that v{si) = v{sf) where Si is the first switching time 
and s f is the time needed to steer the last switching point to 
the south pole. 

Proof of the lemma. Consider the problem of connecting 
the south pole to the north pole in minimum time through 
the system 



F'siz) + uG'siz) 



(20) 



where z e and -Fs(z), G'giz) = -Gs{z). 

The trajectories of system (20) coincide with those of the 
system (4)-(9), but the velocity is reversed. Therefore the 
optimal trajectories for the new problem coincide with the 
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Fig. 6. Comparison between tlie optimal strategy for our system and in 
the RWA 



optimal ones for the system (4)-(9) connecting the north pole 
to the south pole, and the time between two switchings is the 
same. 

Moreover, if we perform the change of coordinates 
{zi,Z2,Z3) (yi, 2/2, 2/3) = {-zi,Z2,-Z3), then the new 
problem becomes exactly the starting problem, and so we 
deduce that, if we have more than one switching, it must be 

V{S^) = v{sf). 

Proof of Proposition 3. First one can easily see that the 
only possible trajectories steering the north to the south pole, 
with only one switching are those described by the proposi- 
tion. So, since by Proposition 2 the total time for trajectories 
with more than two switchings is larger than 27r, it remains to 
compare our candidate with the trajectories containing a sin- 
gular arc and with those with exactly two switchings. In the 
first case the trajectories must be of the type BtSgBt' and the 
only possibility ist~t'= tt— arccos(cot^(Q!)), while the to- 
tal time is 27r— 2 arccos(cot^(a))+2 arccos(cot(a))/ 003(0;) 
which is larger than 2tt. 

One can observe that v{.) is an increasing function if a > 
7r/4 and therefore, if we apply Lemma 1, we obtain that for 
an optimal trajectory with more than one switching it must be 
Si = Sf. In particular the bang-bang trajectories with exactly 
two switchings joining the north pole to the south pole and 
with Si = Sf < IT can be explicitly determined and their 

corresponding total time is 27r + 2 arcsin (^ 2sin(a) ) ■ 

Proof of Propositions 4 and 5. If a < 7r/4 then 
v{0) = v{Tr) = TT, moreover v{.) is increasing between and 
arccos(— tan^ a) and decreasing between arccos(— tan^ a) 
and TT (see Figure 1). 

Therefore, given a < tt/A and s £ [0,7r] with s 7^ 
arccos(— tan^ a), there exists one and only one time s'{s) G 
[0,7r] different from s, such that v{s) = v{s'{s)). 



Notice that s and s'(s) satisfy the nice property 

s + s'{s) ^v{s). (21) 

Indeed both s and s'{s) satisfy the following equation in 

t G [0,7r]: 

/I , A sinm 

cot -v(s) = , , n =^ 

\2 J cos(t) +cot^a 

cos ^iw(s) — <^ = — cos ^2'^'^'^''^ ""^^^ 

Therefore, since |f (s) — t G [— tt, tt] Vs, t G [0, tt] and 
s'(s) 7^ s, it must be: 

s'{s) — 2^'-^'' ~ 2^^'*'' ^ * ^ * ^ ^ ''^i^)- 
So we deduce that there are two possible cases: 

(*) Sf = s'(Si) 
(4) Sf = Sj 

The description of candidate optimal trajectories is simpli- 
fied by the following Lemma, of which we skip the proof. 
Lemma 2: Set: 

, / cot (ii;(s)) — sin(Q;) 

Z{s) = - -cot(iw(s)) 

^ \ sin(a) 



where p= ycot^ (^^^'~*^) +sin^(a) 

Then, if e{s) is defined as in (15), and X+ Fs + Gs, 
:= Fs — Gs , we have: 

^0{s)Z{s) ^ gi'(s)X-gt,(s)X + 

Notice that the matrix Z{s) G so(3) is normalized in 
such a way that the map t t-^ e*-^(*) g 50(3) represents 

a rotation around the axes R{s) = sin(a) with 

V cot(iz;(s)) ; 

angular velocity equal to one. 

Let us study the two possible cases described above: 

(Jit) Suppose that the optimal trajectory starts with u ^ —1 
(the case it = 1 is symmetric) and has an even number n of 
switchings. Then it must be 

S = e^/^" e''(^')^:. . . . e^(^-)^^ e^'^" N (22) 
V ' 

n — l times 

where N and S denote respectively the north and the south 
pole, and we have that 

from which we deduce that Sj must satisfy 

—n6{si) = TT + 2pTT for some integer p. 

It is easy to see that a value of Si which satisfies previous 
equation with p > Q doesn't give rise to an optimal trajec- 
tory (since, roughly speaking, the corresponding number of 



switchings is larger than the number of switchings needed 
to cover the whole sphere). Therefore in previous equation 
it must be p = 0. 

If n is odd the relation (22) becomes 

S = e'f^^ e-(^')^:. . . . e"(^-)^^ e^'^'iV (23) 
« ' 

n—1 times 

and, moreover, by simmetry: 

Then, combining with (23) and using the relation (21), we 
find; 

" V ' 

2n times 

Since e'^^^ N is orthogonal to the rotation axis R{si) 
corresponding to Z{si), previous identity is satisfied if and 
only if nO{si) ~ irrni with m positive integer As in the 
previous case, for an optimal trajectory, it must be m = 0, 
and therefore the proof of Proposition 4 is complete. 

(<|k) For simplicity call si = Sf = s. Assume, as before, 
that the optimal trajectory starts with m = — 1 . If this 
trajectory has n = 2g + 1 switchings then it must be 

In particular the points e^'^^^ S and e**^ N must belong 
to a plane invariant with respect to rotations generated by 
Z{s) and therefore the difference e^-^ N — e~*^ S must 
be orthogonal to the rotation axis R{s). 
Actually it is easy to see that this is true for every value 
s G [0,7r], since both e~^^ S and e^^ N are orthogonal 
to R{s). Moreover, since the circle passing through e'*^ N 
and e"''^ S corresponding to the rotations around R{s) 
has radius 1, it is easy to compute the angle (3{s) between 
these points. In particular the distance between e'''^ N and 
g-sX+^ coincides with 2sin(^^^^) , and so one can easily 
obtain the expression: 

f3{s) = 2 arccos(sin(a) cos(a)(l — cos(s))). 

Then Proposition 5 is proved when n is odd. 

Suppose now that the optimal trajectory has n = 2q + 2 

switchings, then we can assume without loss of generality 

that 

S ^ e'^x- ^v{s)x+ ^qe{s)z{s)^sx- jy ^ 

First of all it is possible to see that e~'"(*'^^e~^^ S is 
orthogonal to R{s). So it remains to compute the angle /3(.s) 
between the point e^-^ N and the point e~"('*)"^ e~^^ S 
on the plane orthogonal to R{s). As before the distance 
between these points coincides with 2sin(^^^). 
Instead of computing directly l3{s) we compute the difference 
between the angles /3{s) and the angle (3{s). We know that 

2sin(M_M) ^ \e~^i^)x\'^^- S - e-^^^5| = 



Using the fact that s and s'{s) satisfy the relation v{s) = 
v{s'{s)) one can easily find that 

\e-'^'S- e^'^'^^^S] = 2^1-cos2(a)sin2 Q^Ks)^ 
Therefore 

(3{s) = f3{s) — 2 arccos ^cos(q;) sin ^^"(*)^ ^ ■ 

This leads to f3{s) — /3{s) = 0{s)/2 and the proposition is 
proved also in the case n is even. 
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